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Abstract 

We present variational Monte Carlo (VMC) results for the Gutzwiller-projected 
coexisting state including both antiferromagnetic (AFM) order and supercon¬ 
ducting (SC) order in the two-dimensional t — t' — t” — J model. By further con¬ 
sidering off-site spin correlation between electrons, in contrast to earlier VMC 
results [Phys. Rev. Lett. 102, 027002 (2009)], we find the apparent competi¬ 
tion between AFM order and SC order near the underdoped regime instead of 
coexistence. The local ferromagnetic correlation introduced by spin-spin Jas- 
trow correlators seem to be responsible for the disappearance of the coexisting 
state. We also demonstrate that the quasiparticle spectral weight from upper 
(lower) AFM band are strongly diminished (enhanced) by the spin-spin corre¬ 
lation. This result obviously leads to the loss of antinodal electron pockets and 
the appearance of nodal hole pockets as passing from the AFM phase to the 
SC phase in hole-doped cuprates, which is in consistent with the observation by 
angle-resolved photoemission spectroscopy. 

Keywords: Variational Monte Carlo method; Strong electron correlation; 
Cuprate superconductor 


1. Introduction 

The doping phase diagram near the underdoped region is one of the impor¬ 
tant and long-debated issues with high-Tc cuprates [T]. Since the parent com¬ 
pound is antiferromagnetic (AFM) Mott insulator, the AFM correlation plays 
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a significant role in the emergence of superconductivity by doping charge car¬ 
ries. The intrinsic proximity of the superconducting (SC) phase with the AFM 
phase is also shared by the phase diagrams of other SC materials, such as iron 
pnictides [2] and heavy fermion superconductors [3]. The multi-layered cuprate 
superconductors exhibit the coexistence of AFM and SC states at underdop¬ 
ing discovered by nuclear magnetic resonance measurements 00. However, 
the AFM phase and the SC phase never coexist in the phase diagram of single¬ 
layered cuprates such as La 2 -a;Sra;Cu 04 |5] and Bi2Sr2Cu06+5 [7]. In particular, 
Bi2(Sr2_a;Laa;)Cu06-i-5 systems shows that the three-dimensional AFM region, 
separated by the SC phase, even survives until a high underdoping level [5]. 

The existence of the coexisting state has been found by analytical and nu¬ 
merical approaches in Hubbard—type models 0 Uni EU [13 [13 m El] and 
t — J—type models [13 EB E3 Ea 123 EH [22] • They seem to contribute the un¬ 
derlying mechanism to the coexisting state observed in multi-layered cuprates. 
However, a proper mechanism to explain why these two phases do not like to 
coexist in single-layered cuprates remains needed. Interestingly, some previous 
studies proposed the spin-bag mechanism for superconductivity since two spin 
bags would attract each other to form a Cooper pair and lower the total energy 
[23 El mi. As for doping more holes, therefore, it is necessary to re-examine 
how the local distortion of the AFM background around holes influences AFM 
order and SC order. 

On the other hand, one of the most exciting experimental results is the 
observation of quantum oscillations in the hole-doped cuprates which pointed 
to electron pockets [23122]. In particular, they proposed that these electron 
pockets probably originate from the Fermi surface reconstruction caused by 
the onset of a density-wave phase, e.g. the AFM phase. Unfortunately the 
electron-like Fermi pockets have never been found in most of hole-doped cuprates 
using angle-resolved photoemission spectroscopy (ARPES) [23 123 133 ■ Thus, 
to comprehend the loss of the electron pocket observed by ARPES experiments, 
we inquire to what extent into the electronic correlations ignored in mean-field 
calculations. 
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In this work, we study Gutzwiller’s trial wave functions with the coexis¬ 
tence of AFM order and SC order by means of variational Monte Carlo (VMC) 
method. To improve the trial state, we further consider the off-site correlations 
between two electrons by applying suitable Jastrow correlators. Surprisingly, 
the long-range AFM order is strongly enhanced due to the local ferromagnetic 
(FM) Jastrow correlation, or precisely local AFM distortion, giving rise to the 
disappearance of the coexisting state near underdoping in the phase diagram. 
Besides, the spin-spin correlation in the non-coexisting state transfers the quasi¬ 
particle spectral weight from the antinodal electron pockets to the lower AFM 
band, and also the nodal hole pockets can remain until superconductivity oc¬ 
curs. Therefore, it is expected that the signal of the electron pockets around 
antinodes cannot be found in many hole-doped compounds by using ARPES. 


2. Theory 


Let us begin by the Hamiltonian on a square lattice of size 16 x 16, 


H — + J ^ 

*0.'^ (hi) ^ 


-n,nj 


( 1 ) 


where the hopping tij = t, t', and t” for sites i and j being the nearest, second- 
nearest, and third-nearest neighbors, respectively. Other notations are stan¬ 
dard. We restrict the electron creation operators to the subspace without 
doubly-occupied sites. In the following, the bare parameters , J)/t in the 
Hamiltonian are set to be in the hole-doped regime: (—0.3,0.15,0.3). In order 
to understand how AFM order and SC order compete in variational phase di¬ 
agram, we choose the mean-field ground state including both AFM order and 
SC order (AFSC) as a starting point, 


/ 

I'i’AFSc) = n 7kt7-kt|0), (2) 

k.s—{a.b} 


where the prime means the product only includes momenta inside the magnetic 
zone boundary (MZB). Note that s represents the quasiparticle coming from the 
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upper AFM band (s = b) or the lower AFM band (s = a). The Bogoliubov’s 
quasiparticle operators are defined as 


iL = Kska - O-WkslkiT- 


(3) 


The coefficients and are the BCS coherence factor of AFM quasiparticles 
corresponding to the s band, 




- ( IH-—- 

2 V 

l-K)^ 


(4) 


where the AFM band dispersion (^k ^k = (^k ± Sk+q) /2. 

Here £k is the normal-state dispersion. Ak(= 2A(coskx — cosky)) is d-wave 
pairing amplitude and m AFM order parameter. The annihilation operators for 
AFM bands, Sko-, are given by 



o-/3k 

Qfk 


Ck(7 

*^k+Qcr 


(5) 


with Q = (tt, tt) and the coefficients 


al = 


Pi 



( 6 ) 


In order to introduce more correlations in the mean-field wave function, we 
first formulate the trial wave function fixing the number of electrons P/Ve with 
on-site Gutzwiller projector Pg{= Ili (1 ~ and charge-charge Jastrow 

correlator {Pj^) 

\^cc) = Pn^PgPj^I'^afsc)- (7) 


More importantly, we also consider the correlation between spins by using spin- 
spin Jastrow correlator (Pf'®), 

I'l'ccss) = (8) 
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The Jastrow correlator is constructed by classical Boltzmann operator, P} = 
, encoding the intersite correlations. For the sake of simplicity, Hi depicting 
charge {i = CC) and spin {i = SS) parts are chosen to be diagonal in real-space 
configuration. The charge-charge Jastrow correlator describes the short- and 
long-range correlations between holes in the lattice system. Thus, 

Hcc = J2v^jninj, (9) 

i<j 

with rjij = Here is the chord length of [ry — rj\ and = 

1 — '^crfiia- We consider three parameters the nearest (7 = 1), second- 
nearest (7 = 2) and third-nearest (7 = 3) neighbors, standing for short-range 
hole-hole repulsion if < 1. The factor r“ denotes attractive long-range 
(rij > 1 ) and repulsive short-range {vij < 1) correlations between holes if a > 0 . 

A similar formalism to the spin-spin correlation has been considered at half¬ 
filling [33]. We further imitate the formalism described above to write down the 
spin-spin Jastrow correlator, 

Hss = ^ Sj , ( 10 ) 

i<j 

where Ky = hi{rfjWj ^and the spin operator along 2 : direction at site 
i. The only difference from the charge counterpart is that the sign of SfSj 
determines the type of magnetic correlations. In other words, it will be the FM 
(AFM) correlation if S'fS'J > 0 (< 0). In addition to the parameter (3 con¬ 
trolling the long-range spin correlations, we consider the other three parameters 
^^ 7 = 1 , 2,3 for the neighboring spin-spin correlations. For example of the FM case, 
the short-range correlation would be suppressed when < 1. On the other 
hand, the factor control the long-range {rij > I) and short-range (r^- < 1) 
correlations. In the long-range case of /3 < 0, for instance, would decrease 
the FM correlation but conversely increase the AFM correlation. 

In addition to the ground state, we also propose a trial wave function for the 
low-lying excitation of the Gutzwiller-projected coexisting state simply gener- 
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ated by Gutzwiller projecting the mean-field excited state 


I^AFSc) = (7D^ I'^AFSc)- 


( 11 ) 


Here we have applied the particle-hole transformation [331 into Eq.(lll to 


avoid the divergence from the nodes of the mean-field wave function. The 
Gutzwiller-projected excited state with both AFM order and SG order fixing to 
Ne — 1 electrons is written as 



( 12 ) 


Hence we can compute the excitation energies Ek{= (d^ktrl-^l'l'kcr)” (^o|-ff|^'o)) 
for either upper (s = b) or lower (s = a) AFM quasiparticles. Furthermore, 
the quasiparticle spectral weight measured from ARPES can be obtained by 
calculating 


(^Ll^L)(^ol^o)' 


(13) 


Some details in the VMG calculation should be noticed. The boundary condi¬ 
tion we use is periodic along both directions. In order to achieve a reasonable 
acceptance ratio, the simulation consists of a combination of one-particle moves 
and two-particle moves. The variational parameters of the Gutzwiller-projected 
coexisting state are optimized by using the stochastic reconfiguration method 
|36j . All physical quantities are evaluated using the optimized parameters. We 
also take a sufficient number of samples (= 2 x 10®) to reduce the statistical 
errors, and keep the sampling interval (~ 40) long enough to ensure statistical 
independence between samples. 

3. Results 

We first consider the trial state with only the charge-charge Jastrow cor¬ 
relator to better demonstrate the variational phase diagram. Then we further 
include the spin-spin Jastrow correlator to see how the phase diagram changes. 
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Figure 1: (a) Variational phase diagram plotted by staggered magnetization Ms (squares) and 
superconducting order parameter Asc (circles). Filled and empty symbols represent l^cc) 
and I’F([ 7 ( 755 ), respectively, (b) The difference of the energy components between l^'ccss) 
and l^cc) a function of hole doping <5 in 16 x 16 lattice. 


Order parameters shown in the phase diagram are determined by the staggered 
magnetization 

Ms = (14) 

i 

and the long-range pair-pair correlation function 

Cpp{R) = ^ E Ac..a'(AUA,+«,a')- (15) 

The creation operator A| ^(= ~ creates a singlet on the bond 

{i,i + a), a = x,y. The factor Xa^a' describes d-wave symmetry: Xa^a' = 1(~1) 
as a = a'{a ^ a'). 

In Figj^a), without the spin-spin Jastrow correlators as indicated by I'l'c'C'), 
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there exists a region showing the coexistence of AFM order and SC order within 
doping 6 < 0.125 in the phase diagram [ 20 l EU [ 22 ], where Ms and Asc{= 
\/Cpp{R > 2)) are finite. Let us turn to the case with both charge-charge and 
spin-spin Jastrow correlators denoted by I'Lccss)- Obviously the coexisting 
region disappears and a clear boundary separating the AFM phase and the 
SC phase shows up at doping 5 = 0.156. Note that near the boundary the 
spin-spin Jastrow correlator can greatly improve the ground-state energies from 
0.3% to 0.7%. From the numerical optimization, we find the spin-spin Jastrow 
correlator can provide a conduit to vary the mean-field AFM order in 
Surprisingly, the optimized spin-spin Jastrow parameters slightly display short- 
range FM correlations in the AFM background (e.g. at <5 = 0.156 the spin-spin 
Jastrow weights for 7 = 1, 2 and 3 would be increased to 1.12, 1.02 and 1.01, 
respectively). The local FM correlation introduced by the Jastrow factors is 
harmful to the mean-field AFM order. To make them balance, it is inevitable to 
largely enhance the AFM background in The surprising competition 

between AFM order and SC order near the phase boundary is mainly due to 
the hugely enhanced AFM order further leading to the diminished SC order. 

To further demonstrate the energy competition, we analyze the difference 
of the energy components in the Hamiltonian between |d>c'css) and l^'cc) 
shown in Figj^b). Our data clearly show that within 0.04 < <5 < 0.2 the spin 
Jastrow correlator helps the trial mean-field state gain much more energy from 
the second-nearest-neighbor hopping term. On the other hand, the competing 
energy primarily comes from the spin-spin superexchange interaction. From 
real-space point of view, holes prefer to move along diagonal direction in strong 
AFM background so that the hopping energy from the second nearest neighbors 
{t') is likely to compete with the superexchange energy (J). 

In momentum space, it is apparent that the t' energy gain would influence 
how the band dispersion evolves from Fermi pocket to Fermi surface as increas¬ 
ing doping. In Figj^a)-(c), the difference of the momentum distribution func¬ 
tion between |4>c'css) and |4'cc') shows how electrons distribute in the band 
structure. At J = 0.125 (Fig[^a)), obviously electrons in the system would 




Figure 2: The difference of the momentum distribution function between liltccss) s-nd l^cc) 
for doping (a) <5 = 0.125, (b) <5 = 0.156 and (c) <5 = 0.188 plotted in the first Brillouin zone, 
(d) The next-nearest-neighbor energy component, —At' cos{kx) cos{ky). The black diamond 
is the half-filled Fermi surface. White (Purple) regions present the positive (negative) values. 
Red lines mean zero. 


prefer to stay around ’’hot spots” rather than living near nodes and antinodes, 
which hole pockets and electron pockets seem to be observed as well. The hot 
spot is defined as the momenta along the MZB that can be connected by (tt, tt) 
momentum scattering. Once doping is increased to 0.156 which is the phase 
boundary (Figj^b)), hole pockets become larger and electron pockets slightly 
shrink. Now that electrons like to circle just outside the electron pockets, they 
attempt to form a large Fermi surface. Indeed, as further increasing doping 
to 0.188 where the long-range AFM order almost disappears (Figj^c)), a clear 
Fermi surface in which electrons cluster together can be seen. So far, we also 
understand the reason why the system gain much energy from f term since the 
hot spots are located right at the purple region shown in Figj^d). 

In Figj^ we compute the spin-spin, hole-hole and pair-pair correlation func¬ 
tions (already shown in Eq.([T5|)) defined as, 

C'cc(R) = ( 16 ) 
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Figure 3: (a) Spin-spin, (b) hole-hole and (c) pair-pair correlation functions for the optimized 
state l^ccss) denoted by red circle (black square) symbols. There are 40 doped 

holes in 16 x 16 lattice ((5 = 0.156). 


Css(R) = (17) 

i 

The doping density we choose to present is 0.156. Figurej^a) illustrates that the 
spin-spin Jastrow correlators indirectly indnce the stronger AFM background 
showing a constant tail in the staggered spin-spin correlation function which 
implies a clear AFM order. Note that The enhancement of the AFM order 
mainly arises from the mean-field wave function Furthermore, we 

find in Figj^b) that the hole-hole correlation function makes no difference even 
if including the spin-spin Jastrow correlators, except that the short-range part 
becomes less staggered. For spin and charge, there is no correlation for their 
long-range behavior. Finally, we can also see in Figj^c) that as considering 
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Pj^ the pair-pair correlation almost vanishes at large distances so that the SC 
properties is not available. 

Next, it would be interesting to examine the low-lying single-particle excita¬ 
tion spectra near the phase boundary. In Figj^ by applying the ansatz (Eq.(|l2|)) 
to the single-particle excitation, we calculate two quasiparticle band dispersions 
(s = a, b) and their corresponding spectral weight for removing one particle 


defined by Eq.(13|. In order to compare with the excitations with/without 
spin-spin Jastrow correlators Pj^, we plot their excitation energy Ek along the 
high symmetric momenta in Figj^a). In the case where the trial state only 
includes the charge-charge Jastrow factors, its optimized mean-field parameters 
A ^ TO. Due to large d-wave BCS pairing contribution, the dispersions thus 
show convex around the antinodes and almost zero gap between the two bands 
at nodes. Especially, the upper AFM band is beneath the lower AFM band 
near the antinodal regions, and hence there is a clear signal of electron pockets 
arising from the upper AFM band shown in FigQc). 

When further considering spin-spin Jastrow correlators, the optimized mean- 
field parameters to ^ A. Such a huge AFM parameter to gives rise to a typical 
AFM band dispersion and opens a AFM gap between these two bands at nodes, 
as indicated by red circles in Figg a). Interestingly, FigQc) shows that near 
antinodes the quasiparticle spectral weight of the upper AFM band disappear 
and transfer to almost entire lower AFM band (see Fig^b)). In particular, a 
clear hole pocket of the lower AFM band centering around Q/2 is also observed 
in FigQb). The Gutzwiller and Jastrow correlators arising from electronic 
correlation firmly influence the low-lying quasiparticle excitation spectra of the 
mean-field state I^'afsc)- Therefore, the loss of the electron pockets due to 
electron correlations provides a route to figure out why electron pockets have 
never been found in most of hole-doped cuprates measured by ARPES. 
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Figure 4: (a) The quasi-particle excitation dispersion for different optimized states (de¬ 
noted in the legend of (c)) along high symmetric momenta at (5 = 0.156. Empty (Filled) 
symbols represent the upper (lower) AFM band and squares (circles) the trial state \^cc) 
(I^CCSs))' Due to much smaller A than m for the trial state lltccss}! simply plot the 
lower AFM band (red circles) below the Fermi level (pink line) except the nodal regions for 
clear demonstration. The quasiparticle spectral weight are obtained from (b) the lower 
AFM band and (c) the upper AFM band. 


4. Conclusions 

Summing up, by using VMC approach we have studied the coexisting state 
with both AFM order and SC order simultaneously underneath the Gutzwiller’s 
projection and Jastrow correlators. We have thereby re-examined the variational 
ground-state phase diagram and found that the AFM phase competes with the 
SC phase as further considering off-site spin correlations. The reasoning for 
the competition is that the mean-field AFM order is considerably enhanced due 
to short-range FM correlation introduced by the spin-spin Jastrow factors, fur¬ 
ther leading to the vanished SC order. As well, we have first investigated the 
Gutzwiller-projected quasiparticle excitations of the coexisting state. Based on 
the Gutzwiller ansatz, passing through the boundary between AFM and SC 
phases, we have observed the loss of electron pockets near antinodes coming 
from the upper AFM band and the occurrence of hole pockets near nodes aris¬ 
ing from the lower AFM band as long as the spin-spin Jastrow correlators are 
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included. Therefore, such a strongly correlated electron system needs to be 
carefully inspected in the explanation for the low-lying quasiparticle excitations 
observed by ARPES experiments. 
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